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Abstract 

If the particles produced in a nuclear collision undergo collective flow, 
the reaction plane can in principle be determined through a global event 
analysis. We show here that collective flow can be identified by evaluating the 
reaction plane independently in two separate rapidity intervals, and studying 
the correlation between the two results. We give an analytical expression 
for the correlation function between the two planes as a function of their 
relative angle. We also discuss how this correlation function is related to the 
anisotropy of the transverse momentum distribution. 
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1 Introduction 



In searching for evidence of the formation of quark-gluon plasma in ultrarelativistic 
nucleus-nucleus collisions, one is led to address the question whether the matter 
produced in these collisions can be considered, at least locally, to be in thermal 
equilibrium. Local equilibrium implies that the matter behaves collectively, which 
may have observable consequences. When collective flow is present, the evolution 
of the system is determined by the pressure gradient, and is therefore strongly 
influenced by geometrical factors (nuclear size, impact parameter...). Effects of 
geometry on global event shapes have made possible the detection and the study of 
collective flow of nuclear matter at colliding energies up to 1 GeV per nucleon||TJ. 
However, few such studies have been undertaken at ultrarelativistic energies, and 
no conclusion has been drawn so far@. 

Most analyses at intermediate energies have been concerned with the determi- 
nation of the flow direction, which is the direction of maximum kinetic energy flow. 
However, the flow angle (angle between collision axis and flow direction) decreases 
with increasing energy and cannot be measured at ultrarelativistic energies; this is 
because longitudinal momenta are much larger than transverse momenta. When 
it can be measured, the flow direction gives an experimental determination of the 
reaction plane, which is the plane spanned by the collision axis and the impact 
parameter. The latter can also be determined independently by measuring the 
transverse momentum transfer between target and projectile regions 0. Although 
this method has not been successful at SPS||, it has been recently argued that it 
should give results with heavy nuclei[|]]. 

Since the flow direction merges with the collision axis at ultrarelativistic ener- 
gies, whether or not there is collective behavior, evidence for collective flow should 
rather be sought for in the transverse directions. Note further that while analy- 
ses carried out at intermediate energies are mostly concerned with nucleons, the 
large number of mesons created at ultrarelativistic energies in the central rapidity 
region offers an opportunity to study collective flow, not only of spectators and/or 
participant nucleons, but also among the produced particles. Starting from these 
observations, we have proposed a new signature in a previous work||. The idea is 
that for peripheral collisions, the region where nucleon-nucleon collisions take place, 
when projected onto the transverse plane, is anisotropic: it has a smaller size in the 
direction of impact parameter than in the perpendicular direction. This causes the 
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matter produced in the central rapidity region to flow preferentially in the direc- 
tion of the impact parameter, which results in a corresponding anisotropy of the 
transverse momentum distribution. This anisotropy should increase with impact 
parameter. 

We propose to study collective flow in the plane orthogonal to the principal flow 
direction (here, the collision axis). This has also been done at lower energies ||, 
where it was found that matter escapes preferentially in the direction orthogonal 
to the reaction plane. This was referred to as the squeeze-out effect. By contrast, 
we predict a larger flow in the reaction plane. However, this is only an apparent 
contradiction since the effects that come into play are very different . Squeeze-out 
results from an interaction between the participant nucleons, which try to escape 
the fireball, and the spectators which block their path in the reaction plane. But 
at ultrarelativistic energies, the time it takes for the nuclei to cross each other is 
so short that the particles produced in the central rapidity region do not "see" 
the spectators. Anisotropy results from the interaction of particles in the central 
rapidity region among themselves. 

The main problem one usually encounters in global analyses is that fluctuations 
occur due to the finite multiplicity, which may hide collective effects||. In this 
article we would like to show how statistical and dynamical effects can be separated 
in the anisotropy analysis. The idea is that in a given event, one can perform 
two measurements of the reaction plane by doing global analysis with two separate 
subsets of the emitted particles (by selecting for instance the particles produced in 
two separate rapidity intervals). Then the difference between the two results gives 
a direct measure of the dispersion due to statistical fluctuations. 

In section 2 we recall the definitions of the variables we use in the global analysis |J. 
In section 3 we give a general discussion of finite multiplicity fluctuations. The main 
results are contained in section 4 where we calculate, under quite general assump- 
tions, the correlation function between the two measurements of the reaction plane. 
We show how to measure the ratio of dynamical to statistical effects directly from 
the experiment. In section 5 we show how this can be used to isolate the dynamical 
part in the anisotropy. The relevance of this work to current and future experiments 
is discussed in section 6. 
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2 Global analysis 



In flow analyses, an event is characterized by means of global observables describing 
its shape. At ultrarelativistic energies, we are only interested in the transverse 
directions. We therefore define the 2x2 transverse sphericity || tensor S 1 - by 

M 

S tj = Y. w ^) u ^ v ) u A v ) (!) 

where (ui(v), U2(vj) is the unit vector parallel to the transverse momentum of the 
v th particle, w(v) is a weight, and the sum runs over all the particles detected in a 
given rapidity interval. If there is no particle identification, w(u) can be chosen to 
be equal to some function of the transverse energy Exiy) deposited by the particle 
in the calorimeter, for instance u>(z/) = E^{i>), with (3 equal to some real constant. 
Collective flow usually results not only in a larger number of particles emitted in 
the flow direction, but also in a higher energy per particle in this direction, so that 
(3 should be taken positive, for instance (3 = 1. Which is the best weight to consider 
should be determined on the basis of a more careful analysis. However, the results 
which are presented here do not depend on such details. 

S^j has three independent components and is therefore fully determined by its 
two eigenvalues fi and fi (we choose f\ > $2) and the angle 9 between the x-axis 
and the eigenvector associated with fi, with —tt/2 < 9 < tt/2. Instead of fi and f2 
we choose the variables £ and a defined by 



£ = trS L =h + f- 
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S 1 - can then be expressed as a function of £, a and 9 : 

„± _ £ ( 1 + a cos 29 a sin 29 \ , , 

"2V asin29 l-acos29j ^ 

a measures the relative difference between the eigenvalues of S 1 - 1 , i.e. the anisotropy 
of the momentum distribution. As we said in the introduction, collective flow results 
in anisotropy for peripheral collisions^]. Hydrodynamical calculations predict that a 
decreases linearly with the multiplicity or the transverse energy (which are measures 
of the impact parameter). The highest value of a, obtained for very peripheral 
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collisions, is about 0.25 to 0.3 for a Pb-Pb collision, and slightly less, about 0.2, for 
a S-W collision. 

There is a larger flow in the direction of the impact parameter than in the 
direction perpendicular to the reaction plane 0. Thus we expect that the principal 
axis associated with the larger eigenvalue fx is the direction of the impact parameter. 
Since S is directly measurable, this in turn gives an experimental measurement of 
the orientation of the reaction plane. However, this estimate is reliable only if the 
measured anisotropy originates from collective flow. Statistical fluctuations related 
to the finite multiplicity also generate anisotropy, which must be disentangled from 
the dynamical anisotropy created by collective behavior. It is the purpose of this 
article to show how this can be done experimentally. 



3 Finite multiplicity fluctuations 

Macroscopically, a collision between two spherical nuclei is fully characterized by 
the colliding energy and the impact parameter. However, particle emission is the 
result of microscopic processes (parton-parton or nucleon-nucleon collisions) which 
are not described by these macroscopic variables. For fixed energy and impact 
parameter, this results in statistical fluctuations of global macroscopic quantities 
such as the sphericity tensor S^j. Now, since particles are created independently at 
different points in the system, it is reasonable to assume that S^ is the sum of a 
large number of independent random contributions. Then the central limit theorem 
states that the probability law of S£j is gaussian. To avoid having to deal with too 
many indices, we rearrange the three independent components of S^ into a 3- vector 
S whose components Si are defined as 

S± = S-^i + S22 — £ 

5 2 = S^- S& = Sa cos 26 (4) 

5 3 = 2S£ = £asin20. 

The most general gaussian probability distribution for Si is of the form 
d P 

dS^ K - p [-(' S -' 5 ) T "(5-S)/2] (5) 
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where S = (S) is the average value of S and T is the 3x3 covariance matrix defined 
by 

Tij = (SiSj) - SiSj. (6) 

If one considers two independent random variables S and S' with gaussian proba- 
bilities such as (0), with respective covariance matrices T and T', the probability of 
S + S' is gaussian, with a covariance matrix equal to T + T'. Thus T^- is proportional 
to the number of particles M used in the analysis. 

Since the system is symmetric with respect to the reaction plane for spherical 
nuclei, so must be the probability (|5|). We choose x to be the direction of the impact 
parameter. Then the symmetry with respect to the reaction plane changes 9 into 
—9. Eq.(0) is invariant under this transformation if 

S 3 = 

T 13 = T 23 = 0. (7) 
Thus S depends only on two parameters 

S x = £ 

S 2 = £a. (8) 

at represent the anisotropy associated with S, that is the anisotropy in the emission 
law. Although S is the average value of S, a is in general not equal to the average 
value of a, (a) : for an isotropic emission, a — 0, but in general a > with 
a finite number of particles. While at represents the anisotropy associated with 
macroscopic effects, for instance with collective flow, the average value (a) also 
receives a contribution from finite multiplicity fluctuations. 

If the emission is weakly anisotropic, i.e. if a 1, the covariance matrix T is ap- 
proximately the same as for an isotropic distribution. For an isotropic distribution, 
the probability @ must be invariant under rotations, that is under transformations 
9 — » 9 + 9 with #o fixed. This implies 

T 12 = 

T22 = T33. (9) 

Then the probability (|^) can be integrated over £ with the result 

dP 4a ( a 2 + a 2 - 2aa cos 2#\ . . 

r exp (10) 



dad# no 2 \ a 
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with 

a=^. (11) 

a is the order of magnitude of the anisotropy created by statistical fluctuations 
alone, as can be seen easily if at = 0: the probability distribution for a, apart from 
the preexponential factor a which arises from the jacobian transforming {a, 8,9) 
into (Si, S2, S3), is a gaussian of width a/y2. Note that since both T 2 2 and 8 are 
proportional to the number of particles M used in the analysis, a scales like 1 / yM. 
For an uncorrelated emission of identical particles, one finds a = (1/ y/M) (w 2 ) / (w) 2 , 
where w is the weight of the particle in the sphericity tensor, Eq. ([!]). We have chosen 
to normalize Eq. (|T0|) to unity in the interval [0, ir/2] (rather than [— 7r/2, 7r/2] because 
the function is even), that is Jq^ 2 d9 f$ da(dP/dad9) = 1. This will be checked below. 

It is convenient to express the probability as a function of the scaled quantities 

X = a/a and x — (12) 

Eq.fliDD becomes then 

= ^exp (-x 2 - X 2 + 2 X X cos 29) (13) 

This equation only involves the dimensionless parameter x, which will play a crucial 
role in our analysis. Physically, x represents the ratio of the anisotropy a generated 
by dynamical collective effects to the typical anisotropy o yielded by statistical 
fluctuations. Note that strictly speaking, Eq.([T3|) also involves the parameter a 
since x varies from to 1/cr. However, d < 1 for a large system and the probability 
( |i3|) decreases exponentially for x ^> 1. Thus we will let x vary from to +00 when 
we integrate over x- 

Eq . (|13|) can be integrated over 9 using the modified Bessel function Jo defined in 
Eq. ( |A.2| ) . One gets then the distribution of the scaled anisotropy x : 

dP 

— = 2 X exp(-x 2 - X 2 )h{2xx)- (14) 
dx 

To check that this distribution is normalized to unity, we integrate (fbl| ) by parts 
and use the relation I\{z) = dl /dz. Then the formula 



exp(- X 2 )Wxx)dx = V exp(x72)/„/2(r/2) (15) 
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with v — 1 gives the result, using the fact that h/2(z) = 2 sinh z/\ / 2itz. 

In principle, the distribution (|14D can be compared to the experimental distri- 
bution of the anisotropy a: \ and the scale factor a can be fitted so as to obtain 
the best agreement with the data. However, this may be difficult in practise, as we 



are going to see shortly Fig.l displays (l/x)(dP/dx) (that is the distribution (|Tj) 
divided by the factor x arising from the jacobian), for three values of the parameter 
X- For x — (no collective flow), this quantity would be a gaussian of width l/y/2 
centered at x = as can be seen from Eq.(|14|). On the other hand, if x > 1; one 
easily shows that (l/x){dP/dx) reaches its maximum at a non vanishing value of 
X, which becomes closer to x as x increases. This can be used|| as a signature 
of collective flow. If x < 1 ; however, the maximum is reached at x = an d the 
distribution is very close to a gaussian as illustrated in Fig.l for x — 0.5. Thus, the 
effect of x is simply to increase the width of the ^-distribution, compared to % = 0. 
Since the a distribution has the same shape as for x = 0, it is impossible to extract 
X from the anisotropy distribution alone. 



Eq.(|l3D can also be integrated over x, which yields the probability distribution 
of the angle 6 

dP 2 

— = - exp(-x 2 ) {l + V^X cos 26 [1 + erf (x cos 26)} exp (x 2 cos 2 26) } (16) 
where 

erf(x) = ^= f e" i2 dt (17) 

V7T JO 

is the standard error function. dP/d6 is a decreasing function of 6. If statistical 
fluctuations are large compared to the dynamical anisotropy, that is if x ^ 1> 
Eq.flT6"D reduces to 

^ = l(l + ^ X cos2d)+0(x 2 )- (18) 

The anisotropy results here in a small deviation, with amplitude proportional to x, 
from the constant value corresponding to isotropic emission. On the other hand, in 
the limit where x ^ 1 (strong anisotropy), Eq.([T6|) becomes 

^i§e X p(-4 S V). (19) 
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In this case, the probability is a gaussian of width l/(2\/2x) < 1, centered at 9 = 0. 
Note that Eq. ( |T6|) is of little practical use : 9 is measured from the x-axis which we 
have chosen, in this section, to be the direction of impact parameter x. But since 
there is no direct access to this direction in the experiment, 9, unlike a, is not an 
observable. We do not get any physical information from 9 unless it is correlated 
with another independent evaluation of the reaction plane. 



4 Reaction plane correlations 

If one measures the reaction plane from S 1 - in two separate rapidity intervals with the 
same multiplicity, one obtains two angles 9\ and 92 (measured from an arbitrary fixed 
direction) which are two measurements of the reaction plane. The determination is 
reliable only if 9\ and 9 2 are strongly correlated. In this section, we calculate the 
probability distribution dP corr /d9 of the relative angle 9 = 9 1 — 9 2 . If this probability 
is flat, 9 1 and 9 2 are uncorrelated and no conclusion can be drawn concerning the 
occurrence of collective flow. We expect this to be the case if anisotropy is small or 
statistical fluctuations are large, that is if x ^ 1- If> 011 the other hand, x > 1, we 
expect dP corr /d9 to be strongly peaked at 9 = 0. 

We shall assume that the rapidity intervals are well separated, so that there is 
no correlation between them, and consider the two corresponding sphericity tensors 
as statistically independent. Then 9\ and #2 are two independent random variables. 
We further assume that macroscopic quantities (fluid velocity, energy density) are 
invariant under Lorentz boosts along the collision axis[|10|| and postpone the discus- 



sion of this point to section 6. Since the sphericity tensor S involves transverse 
coordinates only, it is also boost invariant. Thus a and a are the same for both 
rapidity intervals and 9\ and #2 have the same probability distribution, given by 
Eq.flTEl). The correlation function is then given by 

dP 1 W 2 r n / 2 dP dP 

mJ d9 2 %-(9 1 )%-(9 2 ) 5(9 -9,-9,) 



d9 2 J-k/2 J-n/2 d9 y ' d9 

1 W 2 dP dP 



where dP/d6* is given by Eq.([l6|). The factor 1/2 normalizes dP COII /d9 to unity 
between and ir/2. 
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The integration can be carried out analytically (see Appendix A), which yields 
the result 



d/ '"" = ? J 2 M - 4 r- 



dO 



{-(1 + x 2 ) + X 2 [cos 29 (Jo + L )(x 2 cos 29) 
+ (/ 1 + L 1 )(x 2 cos20)]} (21) 



where Jo and Ji are modified Bessel functions of the first kind and L and Hq are 
modified Struve functions. Eq.fl2T|) reduces to 



^= = - + t cos 20 (22) 

d9 71 



if X <C 1, and to 



dP„ 



d# V 7T 



Xexp(-2x^ 2 ) (23) 



if x ^ 1- These asymptotic forms can be deduced directly from Eqs.(|18f), (|19| ) and 
(p0|). Note that for small x, the correlations (deviations from a flat probability) 
are of order x 2 - Thus, when statistical fluctuations become larger than dynamical 
effects, correlations decrease very quickly. On the other hand, if x ^ 1, dP cor r/d# is 
the convolution of two identical gaussians ([19]), that is a gaussian of width 1/(2%). 
The correlation function given by Eq.([H]) is displayed in Fig.2 for three values of Xi 
together with the approximations (p2|) and (^3|). One sees that these approximations 
are very good for x < 0.5 and x > 2 respectively. 

A measure of the correlation strength is obtained by forming the ratio of the 
number of events with 9 > 45° to the number of events with 9 < 45°. This ratio 
is equal to 1 if there is no correlation between reaction planes and vanishes if they 
are strongly correlated. Integrating Eq.(^) over 9, one obtains a simple analytic 
expression for this ratio (see Appendix A): 

rn/2 Ap 

/ orr d9 

Ne>A$° _ A/4 d9 _ 1 

aw = r /4 dP ^ d 9 ~ 2 ex p^ 2 ) - 1 ' 1 " ' ; 



o d9 



This quantity is displayed in Fig. 3 as a function of x- As expected, it decreases from 
1 to when x goes from to +oo. One observes that this ratio differs significantly 
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from unity already for modest values of x, so that collective effects should be seen 
easily by measuring reaction plane correlations. Any deviation of dP corr /d9 from a 
constant value can be attributed to collective flow (within our hypotheses), which 
makes this signature less ambiguous than that associated with the distribution of 
anisotropy (see Fig.l and the corresponding discussion in section 3). 

The ratio in Eq.([24j) can be measured directly, and from its value one deduces the 



value of x- One may then check whether Eq.(21) reproduces the observed behavior 
of the correlation function. 



5 Relation to anisotropy 

The correlation between reaction planes clearly does not exhaust all the information 
we get from the sphericity tensor analysis. We also have a measurement of the 
anisotropy a. The value of x one gets from the analysis of plane correlations fixes 
the shape of the x distribution, using Eq . (|i~4"l) . In order to get the a distribution, 
the scale factor a is required (see Eq.fll^D). As we shall see shortly, this quantity 
can be determined through the measured average value of a, which we denote by 
(a). It is directly proportional to o. 

When we studied reaction plane correlations in section 4, only half of the de- 
tected particles (at most) could be used to construct the sphericity tensor since 
we needed two independent evaluations of the reaction plane in each event. On 
the other hand, when measuring the anisotropy distribution, it is better to use all 
the particles detected in the central rapidity region, in order to minimize statistical 
fluctuations. Then the value of x which has been determined from reaction plane 
correlation cannot be used directly in analyzing the a-distribution. Since a scales 
like 1/a/M and a is independent of M, x must be scaled like vM. If, for instance, 
the set of measured particles is divided into two approximately equal subsets for the 
measurement of reaction plane correlations, and then used as a whole for measuring 
the anisotropy distribution, the value of x must be multiplied by \^2. 

Let us now calculate the average value of a: 

{a) = = " sW/Wx (25) 

The denominator of this expression is equal to unity since dPj dx given by Eq. fll^D is 
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normalized to unity. To calculate the numerator, we integrate Eq.(|l4"D by parts and 
use the relation dl\/dz = (Io(z) + 12(2))/ 2 and Eq . (|i~5|) to calculate the remaining 
integrals. The result is 



(X) = \ [(1 + X 2 )/o(x72) + X Z h(rm\ exp(-x72). (26) 

Thus, if one measures \ an d (ct), the last two equations gives the scale factor a and 
thus a = 0% From Eq.([25l) one gets immediately 

M = W (27) 

This quantity is displayed in Fig.4. For small anisotropy, x <S 1, the \ distribution 
( 14|) is gaussian and one gets 

A ^ -t=* (28) 

(a) V 71 " 

In this case most of the observed anisotropy results from fluctuations and the dy- 
namical component a is only a small fraction of the average anisotropy. On the 
other hand, for \ ^> 1, statistical fluctuations become negligible and the average 
anisotropy (a) is close to a. Asymptotically, 



a 1 

(a) ~ 4x : 



- 1 " (29) 



Once a and a are determined, the measured a distribution can be compared to the 
theoretical prediction, Eq . (|H]) . 



6 Discussion 



Let us recall and discuss the hypotheses on which our calculations are based. The 
first hypothesis was made at the beginning of section 3, where we assumed that the 
sphericity tensor Sy has a gaussian distribution. This is true if it can be considered 
as the sum of a large number of independent sources, which is a reasonable assump- 
tion if the nucleon-nucleon collisions creating the particles are incoherent. However, 
deviations from this behavior can occur due to jets, which result in strongly corre- 
lated, strongly anisotropic emission of particles. At very high energies, for instance 
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at LHC, one expects a large number of jets per event, and pairs of jets can be 
considered independent so that our statement holds : the only consequence is that 
the number of independent sources is the number of jets rather than the number of 
produced particles, so that we expect larger statistical fluctuations. On the other 
hand, if only a few jets (say, one or two) are produced in each event, which may be 
the case at lower energies, significant deviations from the gaussian distribution may 
occur. These would cause deviations of the anisotropy and angle distributions from 
the shapes predicted by our model, given by Eqs.([L4]) and (pTQ . 

The second hypothesis was made in section 4 where we treated the sphericity 
tensors measured in two separate rapidity intervals as independent variables. This 
is not strictly true if the rapidity intervals are too close to each other : for instance, 
a resonance decay or a pair of jets can give contributions to both rapidity intervals. 
This can be avoided by taking two remote rapidity intervals. This will be possible at 
RHIC and LHC if detectors have a large rapidity acceptance. In current experiments 
at CERN and Brookhaven, the rapidity window is not so large and one may be 
constrained to work with adjacent rapidity intervals. One may check directly, as a 
test of statistical independence, that the value of x deduced from the analysis indeed 
scales like a/M, with M the number of particles used in the analysis. 



We also assumed in section 4 that the probability distribution of S ± was the same 
for both rapidity intervals as a consequence of Bjorken's scenario for the fluid 



evolution. This scenario is known to be unrealistic in current experiments since 
the measured rapidity distributions are not flat. However, this hypothesis is not 
crucial here. It simplified the calculations since we used the same value of the scaled 
anisotropy x f° r both rapidity intervals. Without this hypothesis, we wouldn't have 
obtained analytical expressions, but the qualitative ideas would remain the same. 
In particular, the correlation between reaction planes could still be used to identify 
collective flow. Note further that recourse to Bjorken's scenario can be avoided for 
a symmetric collision: if the two rapidity intervals are chosen symmetric of each 
other in the center of mass frame, they are equivalent by symmetry and the value 
of x is the same for both. Even more generally, it is reasonable to assume that the 
anisotropy a depends only weakly on the rapidity since it is not much affected by the 
longitudinal expansion ||. If the two rapidity intervals have the same multiplicity, 
the statistical fluctuations a should also be comparable and therefore there is no 
reason for x to change drastically from one rapidity interval to another. 

Although one must be careful in applying our results according to the above dis- 



13 



cussion, reaction plane correlations appear to provide a powerful tool for identifying 



collective flow. The ratio defined in Eq.(24) is very sensitive to collective effects: 



for x — 0.3, this ratio is about 0.85, and such a deviation from unity would be seen 
clearly in an experiment. We have shown |J that the anisotropy distribution alone 
should allow to identify collective flow at AGS and SPS with heavy nuclei (Au or Pb 
projectiles). With the plane correlation method presented here, the possibility is not 
excluded that collective effects could be seen with lighter projectiles, for instance 
with 32 S or 40 Ca beams on heavy targets. 
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Appendix A: Derivation of the correlation function 

Inserting Eq. (|T3"|) in Eq. one may write the correlation function in the form 

g r+oo r+oo . . 

— exp(-2x 2 ) / XidXi X2^X2exp (-Xi - X2) (A- 1 ) 

7T JO JO v 7 

x I"' 2 dBi exp [2xxi cos 26 l + 2xX2 cos 2{6 1 - 6)] 

J-tt/2 

This can be integrated over 6\ by using the modified Bessel function I Q : 

d9i exp (A cos 2Q X + B sin 26 A = ttI (VA 2 + B 2 ) (A.2) 

J-w/2 V ' 

This gives 

r+oo r+oo . . 

exp(-2x 2 ) ^ Xldxi J ^2^X2 exp [-Xi - X2) 



A P 8 r+oo 
LiJ corr / 1\ I 

d6 



dP 

corr 



d^ 7T 



x /o 2 X Jxi + X 2 2 + 2 X iX2 cos 26 (A.3) 
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Introducing polar coordinates in the (xuXz) plane, defined by xi — r cos0, X2 
r sin0, with r > and < <fi < tt/2, the integral over r is of the type 



drr 3 exp(—r 2 )I (ar) 



as may be checked by expanding I . Eq.(|A.3|) thus becomes 




(A.4) 



d9 



7T 

X 



tt/2 



sin . 



1 +x 2 (l + sin 20 cos 29) 



exp (y 2 sin 2<p cos 29 



(A.5) 



Note that the integrand is invariant under cf) — > tt/2 — <ft, which reflects the fact 
that Xi an d X2 play symmetric roles in Eq.( |A.l| ). The integration range can then 
be restricted to the interval [0, 7r/4]. Making the change of variables tt/2 — 2<p —>■ <f>, 
the integral can be expressed in terms of the modified Bessel functions 



h{z) 
h(z) 



2 r/ 2 



TT JO 

2z 



cosh(^ cos 



2z / >7r / 2 

— / sin 2 0cosh(2;cos 
7T Jo 



(A.6) 



and the modified Struve functions[ll| Lo(z) and Li(z) which have expressions similar 
to Iq and I\, with cosh replaced by sinh. After an integration by parts, Eq. ([A.5 ) 
yields the result 

2 



d9 
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;i + r) + r cos 29 (i + cos 29) 



+ (/ 1 + L 1 )(x 2 cos 29)]} 



(A.7) 



identical to Eq.(pT|). Let us check that this expression is normalized to unity when 
integrated from to tt/2. the functions cos26 l /o(x 2 cos 2$) and ii (x 2 cos 26 1 ) are odd 
in COS26 1 and thus do not contribute to the integral. Using the relations I\{z) = 
dl /dz and Lq(z) = dh /dz + 2/n and the definitions (|A.6|) , the integrals of Struve 
functions can be expressed in terms of integrals of Bessel functions, which can be 
found in the literature or calculated directly by power series expansion: 



/•tt/2 

/ d9 Li (x 2 cos i 
Jo 



tt/2 



'COS0/ O (X cos 



-1 + 



sinh x 



w2 



X 2 



(A.8) 
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and 

fTr/2 _ fir/2 



I d9 cos 6 L (x 2 cos 9) = / d<p cos <pli (% 2 cos 4>) 
Jo Jo 

cosh x 2 — 1 



X 2 



(A.9) 



Normalization of the probability ( A. 7 ) follows immediately. 

Let us finally calculate the ratio defined in Eg. (|2"4|) . Since the correlation function 



( |A.7[ ) is normalized to unity, we only need to calculate f^j 2 (dP COII /d6)d6. On this 
interval, all the terms in Eq. (|A.7| ) give a non vanishing contribution. However, using 
Eqs.([A.8[) and (|A.9|) one obtains 



d9 cos2#/ (x 2 cos2#) = -- - d6L 1 (x 2 cos26) (A.10) 

Jtv/A 2 Jir/4 

r' 2 d9 h{f cos 29) = - I"' 2 d9 cos 20 L^x 2 cos 20) (A. 11) 

Jn/4 Jn/4 



nr/2 1 i-k/2 

d6 cos28 1 {x 2 cos26) = / 

/tt/4 2 Jn/4 

and 

Mr/2 i-x/2 

d9I 1 (x 2 cos29) = - / 

t/4 

The terms involving Bessel functions and Struve functions in Eq. (|2i|) thus cancel 
pairwise in the integration and one gets 

fir/2 a p i 

I ^d9= V -eM~X% (A.12) 
Jtt/4 d9 2 

from which Eq.(|24]) follows immediately. 
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Figure captions 



Fig. 1 : The solid lines display the values of (l/x)(dP/dx), given by Eq.flllTj, 
as a function of x f° r three values of the parameter x- The dashed line is an 
approximation to the curve x — 0.5 by a gaussian corresponding to the same average 
value (x) of Xi calculated from Eq.(^fJ) (see section 5). All curves are normalized to 
unity: f +oo (l/x)(dP/d X )dx = 1- 

Fig. 2 : Solid lines: correlation function defined by Eq.(|2~l~D for three values of \. 
Dashed line: Small x approximation, Eq.(^). Dot-dashed line: large x approxima- 
tion, Eq.pl). 



Fig. 3 : Ratio defined in Eq.(24) as a function of the dimensionless parameter x- 

Fig. 4 : Solid line: ratio of the "dynamical" anisotropy a to the measured average 
value of a, denoted by (a), as a function of % Dashed line: small x approximation, 
Eq . (p8|) . Dot-dashed line: large x approximation, Eq . (|29|) . 
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